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Abstract 

The partial widths are determined for stop decays to top quarks and gluinos, 
and gluino decays to stop particles and top quarks (depending on the masses 
of the particles involved). The widths are calculated including one-loop SUSY- 
QCD corrections. The radiative corrections for these strong-interaction decays 
are compared with the SUSY-QCD corrections for electroweak stop decays to 
quarks and neutralinos/charginos and top-quark decays to stops and neutrali- 
nos. 
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1 Introduction 



The top and stop particles form a complex system in supersymmetric theories. The strong 
Yukawa coupling between top/stop and Higgs fields gives rise to a large mixing of the L 
and R stop states II and £r, which are associated with the left and right chiral top-quark 
states. The mass splitting between the stop mass eigenstates t\ and £2 can therefore be 
quite large. In fact, it is possible that the mass mi of the lightest stop state is even 
smaller than the top mass m t itself [|TJ. 

Depending on the mass values of the particles involved, quite different decay scenarios 
will be realized in the stop-top sector. If stop particles are very heavy, they can decay 
into top quarks and gluinos, 



tj->t + g 



mi. > m t + m g 



In this paper we generalize the analysis of Ref. 0, in which the decay widths for the 
squarks related to the light quark species (q — u, ...,b) were calculated, to the stop-decay 
processes ([I]) in next-to- leading order SUSY-QCD; in this more complex case the stop 
mixing and the non-zero top-quark mass in the final state must be taken into account. In 
a similar way we analyze the crossed channel 



tj and c.c. 



mg > m t + mi. 



(2) 



in leading and next-to- leading order. For small stop masses, in particular for ti, the decay 
channel (|1|) is presumably shut kinematically and decays to quarks and light neutralinos 
or charginos (ti — > txi,bxt) will be dominant |3|].[] For the sake of comparison with 
Refs.p, |^], we re-analyze also these decay modes in next-to-leading order SUSY-QCD. 
Finally, in the exceptional case m t > m t - , the interesting decay mode t — > tixl ma Y 
occur, see e.g. Ref. 0; the partial width of this non-standard top decay has recently been 
determined in next-to-leading order SUSY-QCD in Ref.j7j. 



2 Theoretical Set-up 

To lowest order the partial widths for the stop and gluino decays ([I]) and (0) are given 



r(fi, a ->ts) 



2a s K 



3m 



m? 



ft 



1,2 



h,2 
8m? 



E l,2 



m. 



m 



~ ± 2m t m g sin(2#) 



rn 



*1,2 



— m\ — m~ ± 2mtmg sin(2#) 



(3) 
(4) 



1 For higher-order electroweak stop decays we refer to the recent paper Ref.Q. 

2 As usual, k = (YJj mj — rnfrnj) 1 / 2 , the sums running over all particles involved in the decay 

process. 
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Here m% and are the eigenvalues of the stop mass matrixP] |], ||: 



M 2 RL M 2 RR 



m % + m t + (| _ I s ™) m z cos(2/3) -m t (A t + /i cot (3) 

—m t (A t + fi cot (3) mjj + + \s 2 m m 2 z cos(2/3) 



(5) 



The quantities mQ,mu, /z, and A 4 are the usual soft SUSY-breaking mass and trilinear 
parameters, mz and s w are the Z-boson mass and the weak mixing angle, and tan f3 is the 
ratio of the two vacuum expectation values in the Higgs sector. The diagonal entries of the 
stop mass matrix correspond to the L and R squark-mass terms, the off-diagonal entries 
are due to chirality-flip Yukawa interactions. The chiral states and t R are rotated into 
the mass eigenstates ho and t 20 



ho \ _ ( cos^o sin^o \( t L 
ho J V — sin^o cosfi'o / V 



(6) 



by these Yukawa interactions. The mass eigenvalues and the rotation angle can be calcu- 
lated from the mass matrix (El): 



2 2 1 

*10' *20 2 



M 2 LL + M 2 RR =F [(M 2 LL - M 2 RR ) 2 + 4(M 2 LR ) 2 } 1/2 1 (7) 



sin(2g )= 2 2MIr 2 and oos(2flb) = ^ M ** (8) 

ho *20 *10 *20 

By definition we take ho to correspond to the lightest stop state. 

The mixing angle is an observable quantity, as evident from the decay widths (3) and 
(4). Since the widths of supersymmetric particles are notoriously difficult to measure, 
production processes may instead be adopted for the operational definition of the mixing 
angle in practice. Pair production in e + e~ collisions, e + e~ — > trfj, lends itself as a 
convenient observable |]TT] . 

SUSY-QCD corrections, as exemplified by the diagrams of Fig.|l|(a), modify the stop 
mass matrix and the fields, necessitating the renormahzation of the masses m~ = m~ + 

Sm~ and of the wave functions ho — Zij 2 tj] the renormahzation of the mixing angle can 

be related to the renormahzation matrix Z 1 ! 2 . Without loss of generality, we may assume 
Z}( 2 = —Zo! 2 since, as will be shown later, the reduced self-energy matrix Ej 7 7 (m 2 ~ —m 2 ~ ) 

is antisymmetric. In that case the renormahzation matrix Z 1 ! 2 can be written to order a s 

1 /2 

as the product of a diagonal matrix with the elements Z-- = 1 + 5Zjj/2 and a rotation 
matrix parametrized by the small angle 58: 

1/2 f l + 6Z n /2 \f cosSd sm59\ 

1 1 + SZ 22 /2 J \ -sin 59 cosSOj { > 



The sign conventions follow the SPYTHIA program ||, which is based on Ref.[]lO[ 
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with 86 = Z\ 2 = —Zl[ 2 . Within this formalism the rotational part of Z 1 ! 2 can be 
interpreted as a shift in the mixing angle, 6q i— > 6q—86 = The remaining renormalization 
constants are fixed by imposing the following two conditions: 

(i) The real part of the diagonal elements in the stop propagator matrix Aij(p 2 ) develops 
poles for p 2 —>■ mf , with the residues being unity. This requirement fixes the counterterms 

8m 2 and the diagonal wave-function renormalization 8Z~~. The so-defined renormalized 

stop masses mf are called pole masses. 

l j 

(ii) We define a running mixing angle 0(Q 2 ) by requiring that the real part of the off- 
diagonal elements of the propagator matrix Aij(p 2 ) vanishes for a specific value p 2 = Q 2 
of the four-momentum squared. In this scheme the [real or virtual] particles t\ and £ 2 
propagate independently of each other at four-momentum squared Q 2 and do not oscillate. 

The dependence of the renormalized mixing angle 6 on the renormalization scale Q is 
indicated by the notation 8(Q 2 ). Different choices for Q 2 are connected by a finite shift 
in 6(Q 2 ), which is calculated in the next paragraph. Quite often the renormalized mixing 
angle is defined in an 'on-shell renormalization scheme' |5|, |T2] in which the renormalization 
of the mixing angle cannot be linked to the stop wave functions any more. Even though the 
definitions of the mixing angle are different in the two schemes, the physical observables, 
widths and cross-sections are of course equivalent to 0(at s ). 

In carrying out this renormalization program we find the following expressions for the 
various counter terms: 



8m\. = Re %(m|) 8Z j3 = -Re %(mf.) 86 = -Re S 12 (Q 2 )/[m| - m|] 



Here E^Qo 2 ) and £jj(p 2 ) = dY*ij(p 2 )/dp 2 denote the unrenormalized self-energy matrix 
and its derivative [see also Ref. |12|, ||] : 

Si 2 (p 2 ) = -27rCV« s js 4e ~ A(m t ~ 2 ) - A(m~ tl ) + %m~ g m t c 2 Q B(p 2 ,mg,m t )^ (10) 
S 21 (p 2 ) = E 12 (p 2 ) (11) 
EnGo 2 ) = -4nC F a s { (l + c 2 2§ ) A{m h ) + s^K 2 ) 

- 2 A(m~ g ) - 2 A(mt) ~ 2 (f + m?J B{p 2 , A, m h ) 

+ 2 (p 2 - m~ - m 2 t + 2mgm t s 2 g) B(p 2 ,m g ,m t )^ (12) 
S 22 (p 2 ) = -47rC F a s | (1 + c 2 2§ ) A(m h ) + s 2 2§ A(m h ) 

- 2 A(m~ g ) - 2 A{mt) - 2 (p 2 + m|) B{p 2 , A, m h ) 

+ 2 (p 2 - m\ - m\ - 2mgm t s 2 g) B(p 2 , m g , m t ) | (13) 
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[We have used the standard notation sg = sin# etc.] The first two terms in £12 (p 2 ), 
involving the p 2 -independent 1-point function A, follow from the third Feynman diagram 
in Fig.^(a), the remaining term given by the 2-point function B, corresponds to the 
second diagram in Fig.|l|(a). As both A and B are ultraviolet (UV) divergent^, also 59 is 
UV divergent; in n = 4 — 2e dimensions: 



59 



C F a s 



div 



$7r(m~ — m~ ) 



M0 



J e 



(14) 



The change of the renormalized angle 9(Q 2 ) between two different values of Q 2 is finite: 



9(Q ,2 )-9(Q 2 



167iCFa s mgm t cos(26 l ) 



2 2 
m~ — m~ 



Re 



B(Q' 2 , nig, m t ) - B(Q 2 , m- g , m t ) 



(15) 



This shift is independent of the regularization scheme. For illustration the normalized 
shift relative to 9(m 2 ^): [9(Q 2 ) — 9(m 2 ^)]/9(m 2 ^), is shown in Fig.g. 



3 Stop and Gluino Decays 

The diagrams relevant for stop and gluino decays are presented in Fig.[|]. This set is 
complemented by the self-energy diagrams for the stop particles, the gluinos, and the 
top quarks, displayed in Fig.Q. The Born diagrams are presented in Fig.|^(a) for the 
two decay channels, the vertex corrections in Fig.|3](b), and the hard-gluon radiation in 
FigfJ(c). In Figs.0(b) and (c) only the stop-decay diagrams are depicted, since gluino- and 
stop-decay diagrams are related by crossing. The ultraviolet divergences are regularized 
in n dimensions, infrared and collinear divergences^ by introducing a small gluon mass A. 
The renormalization of the strong coupling constants g s and g s are carried out in the MS 
renormalization scheme at the charge-renormalization scale hr\ a finite shift |HJ between 
the bare Yukawa coupling g s and the bare gauge coupling g s restores supersymmetry at 
the one-loop level in the MS scheme [see Ref.p[ for explicit verification]: 



9s = 9s 



in V3 



(16) 



[Similarly the electroweak t{tXj and tfixj couplings may be written as ae + bY t + cY b , 
with e = e [1 — a s Cp/(87r)] and Y q — Y q [1 — 3a s Ci?/(87r)] in terms of the electromag- 
netic coupling e and the quark-Higgs Yukawa coupling Y q oc em q ] The heavy particles 
(top quarks, squarks, gluinos) are removed from the evolution of a s (fi 2 R ), decoupled 

4 Thc definitions of the scalar functions A and B can be found in the Appendix. 
Since no ggg three-gluon vertices are involved in the calculation, infrared singularities can be reg- 
ularized by a non-zero gluon mass. This method can also be applied to SUSY-QCD processes including 
ggg vertices of massive gluinos. 
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smoothly for momenta smaller than their masses. The masses of these heavy particles 
are defined as pole masses. In the numerical analyses we have inserted the mass of the 
decaying particle for the renormalization scale [Ir. 

The detailed analytic results for stop and gluino decays are presented in the Appendix. 
In this section we illustrate the characteristic features of the results by numerical evalu- 
ation of a few typical examples, which properly reflect the size of the SUSY-QCD effects 
in general. 

The masses and mixing parameters chosen in the examples are calculated from the 
universal SUGRA parameters ||14]| : the common scalar mass m , the common gaugino 



mass m 1( / 2 , the trilinear coupling A , the ratio of the vacuum expectation values of the 
Higgs fields tan/3, and the sign of the higgsino mass parameter fi. The top-quark mass 
is set to m t = 175 GeV, and the Higgs parameter tan/? is fixed to 1.75. From this set, 
the pole masses of the charginos, neutralinos, gluinos, and squarks, as well as the squark 
mixing matrices can be calculated. We use the approximate formulae implemented in 
SPYTHIA ||. We define a degenerate squark mass by averaging over the five non-top 
flavors. The renormalized stop mixing angle 6(m 2 ^) is defined in this basis by imposing 

the lowest-order relation Eq. (8) in terms of the renormalized (pole) masses.0 The mixing 
angles at other renormalization scales can be obtained by adding the appropriate finite 
shifts [see Eq.flTSD]. 

(a) t 2 t + g 

for m = 800 GeV, A = 200 GeV, // > 
In Fig.f|(a) the stop and gluino masses are given as a function of the common gaugino mass 
7711/2. F° r the indicated set of parameters the decay t 2 - ► tg is the only strong- interaction 
decay mode that is kinematically allowed. 

In Fig.f|(b) the width of t 2 is presented in lowest order and in next-to-leading order as a 
function of . Since in this example rrig rises faster than with increasing m^, the 
width drops to zero when mi is increased. The radiative SUSY-QCD corrections vary 
between +35%, at the lower end of the spectrum, and ~ +100% at the upper end of the 
spectrum, i.e., the corrections are large and positive. 

(b) g tj+ttf 

for m = 400 GeV, A = 200 GeV, // > 

In a form analogous to the preceding example, masses and widths are displayed in Figs.[|(a) 
and (b). Since m~ g rises faster than mj , the width increases with increasing gluino mass. 
However, in contrast to stop decays, the SUSY-QCD corrections to gluino decays are only 
modest and negative [~ —10%]. This result is familiar from the analysis in Ref. ||, where 
it had been demonstrated analytically that negative it 2 terms, arising from the crossing 

6 Equivalcntly the mass-matrix parameters tuq, mjj, and A t may be chosen as basic parameters, and 
0(m~ ) and m ( - may subsequently be defined in this basis by imposing the lowest-order relations Eqs. (7) 
and (8). We have checked that the two schemes are equivalent. 
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of diagrams, give rise to destructive interference effects such that the overall correction is 
small. 

(cjt^t + Xj b' = l,-] and h-^b + xt [j = l,2] 

for m = 50 GeV, A = 100 GeV, n < 

For the sake of comparison we have re-analyzed stop decays into neutralinos and charginos. 
The results agree with the analysis in Ref.0, and also with the parallel calculation in 
Ref. ||, with which detailed point-by-point comparisons have been performed. For the set 
of parameters chosen in the figures, only two neutralino decay modes are kinematically al- 
lowed. The SUSY-QCD corrections to the decays into neutralinos are small [Fig.||(b)], typ- 
ically less than 10%. The picture is quite similar for the decays into charginos [Fig.^(b)], 
though the corrections are slightly larger as a result of the massless b quarks in the final 
state. 

(d)t — ti+x; b' = i,-] 

for m = 250 GeV, A = 800 GeV, fi > 

As in the previous case (c), this last example is shown merely as a cross-check with 
Ref.0. As shown in Fig.|8](b), the corrections are small. If this top-decay mode is realized 
in nature, the branching ratio for decays into Xi can i n principle be of the order of 4%. 
[Note that T(t -> bW + ) ~ 1.4 GeV for m t = 175 GeV pf]. 

Since the tj/g/t SUSY decay modes could only be illustrated for a specific set of 
parameters, a general programF] has been constructed for generating all relevant decay 
widths in the stop-top sector. 

4 Summary 

In this paper we have analyzed the SUSY-QCD corrections for stop decays to top quarks 
and gluinos, and gluino decays to stop particles and top quarks. In contrast to earlier 
analyses, the non-zero top-quark mass must be taken into account in these decay modes. 
Moreover, the L/R squark mixing plays an important role. 

We have set up a scheme in which the mixing angle 0(Q 2 ) is defined in such a way 
that the virtual/real stop particles t\ and t 2 do not oscillate for a specific value of the 
four- momentum squared Q 2 . Convenient choices for Q 2 are Q 2 = m~ or m| , depending 
on the problem treated in the analysis. Different conventions are connected by simple 
relations between the associated mixing angles. 

As observed earlier for squarks related to the light quark species, the SUSY-QCD 
corrections are large and positive for stop decays to top quarks and gluinos. They are 
modest and negative for gluino decays to stop particles and top quarks. We have compared 

7 The FORTRAN program may be obtained from plchn@desy.de. 
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these modes with other stop-decay modes, which have been analyzed earlier || ||. In 
contrast to the strong-interaction decays, the SUSY-QCD corrections for electroweak stop 
decays into neutralinos and charginos are small, and so are the corrections for top-quark 
decays to stop particles and the lightest neutralino. 
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6 Appendix 

In this appendix we present the explicit formulae for the calculation of the decay width of t% 
particles to top quarks and gluinos in next-to-leading order SUSY-QCD. The results for the 
decay of ti particles can be derived by interchanging the stop masses and switching the sign 
in front of sin(2#) and cos(2#). The gluino decay width T(g — > itj) can be obtained from 
F(tj — > tg) by adding a factor — m| /(4C^m|) and subsequent analytical continuation from the 

region m^ > mt + m g to the region m g > mt + mr . . The multiplicative factor reflects the 

difference in phase-space, in color/spin averaging, and in the sign of the gluino momentum inside 
the spinor sum [see also Eqs. (3) and (4)]. 

The decay width in next-to-leading order may be split into the following components: 

r NLO = r LO + Re [Ar t + at 9 + AT n + AT V + ad + Ar c + AT f + Ar dec ] 

To allow for more compact expressions we first define a few short-hand notations: 

Ma6c = ml + ml -m 2 c a 2§ = m t m g sm(2<9) 

M = V- 

h 

1/2 



4 i 4 i 4 n/2 2 i 2 2 i 2 2\ 

mr + mx + m+ — 2(mr m~, + m x m, + m^m^ 
q y z v q y q " ' 



where abc = g,t,j with j representing tj. 

We list the components defined above for the next-to-leading order decay width [the 1/e poles 
in the scalar integrals cancel against each other in the final sum]: 

lowest-order decay width: 

IYo = 8N c C F 7ra s (-fx m + 2a 2§ ) N = M \M B \ 2 



8 



top self-energy contribution: 



AT t = ^%^ttq s \M B \ 2 {2(1 - e)A(m t ) + 2A(m g ) - A(m f ) - A{ mi ) 

„2 _ _ \ 1 n/„2 



16iV c C 2 /y 



+ m§B(p t , mj, ) + Ht2§B(p t , m § , m t ~ 2 ) 
- 4m ^ 2 e [ B (Pt > m s . m t~! ) - # (p? , m 5 , m t - 2 )] 

+ 2m 2 [ mi B(pl m- g ,m h ) + ^ gt2 B{pl m s , mtf - Am 2 B(p 2 , A, m t )] } 
vr 2 a 2 c^/x 5t i {A(m~ t2 ) - Aim^) - H g t\B{p\ ,m g ,m~ ti ) + n gt 2B(p 2 , m g , m t ~ 2 )} 



gluino self-energy contribution (for n/ = 6 quark flavors): 

AI> = ^7ra s |.Mb| 2 (n/ - 1){ - A(m g ) + (m\ + m?)B(p?, m 9 ~, 0) + 2m|(m| - m%B{jp% m 9 ~,0)} 



47V" 

|X B | 2 |2A(m t ) - A{m ii ) - A(m~ t2 ) + ^ lgt B(p 2 g ,m- ti ,m t ) + H2gtB{p~,m i2 ,m t ) 
~ 4rn | a 20 |£(p|> , m t ) - B(p 2 g , m~ H , mt)] 



H 7T Tra s 



+ 2mJ [^l-Bfj? 2 ,, m ?1 , m t ) + n~gtiB($\ m~ H , m t )] | 
+ ™s \ M b\ 2 N c {(1 - - 4m|5(pf , A, m 5 )} 
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diagonal stop self-energy: 



Arn =8C F N7ra s \M B \ 2 ^B(p 2 ^,m g ,m t ) - B(p 2 ^,X,m ii ) + 2a 2 gB(p 2 ^,m g ,m t ) 

- H g t\B{p\^m g ,m t ) - 2m 2 h B(p 2 h ,X,m h )} 



[The off-diagonal mixing contribution 



128N C CW 2 2 2 t . .. , 4m 2 m? 

Ar i2 = 2 _ 9 * «s ^ | A ( m * 2 ) " A( * m h > + ~Z B ^h ' m » ' mt > 1 

"\ m t 2 20 

is absorbed into the renormalization of the mixing angle for 9 = 6{m 2 ).] 



vertex corrections: 



AI\, = 64J^TT 2 a 2 N c C 2 F F[ + a 2§ F 2 F + <r 2 ~if 



+ 327V7r 2 a 2 iV 2 C F 
+ 8AAvra s \M B \ 2 F B 



l 26 i 2 -T- u 26 r Z 

Ff + a 2 ~ e Ff + a 2 2§ F£ 



9 



with: 

?F 



Ff = 2(mf + ml)B(p 2 h , m g , m t ) + (m^ + mf + mf)B{p\^ , A, ) 

+ 2(m| - m 2 ^)B{p 2 u A, m t ) - 2m 2 B(p 2 , m g , m t ~ 2 ) - \ra 2 g B{p\m u m iA ) 
+ 4m? (m^ - mfjCip^ , -p t , m t ,m g , ) + 2m 2 (m 2 ^ + m 2 ^ - 2m 2 t )C{p il , -p t , m t , m g ,mz 2 
Ff = - 2B(p 2 h , A, m h ) - 2B(p 2 t , A, m t ) - 45^ , m 5 , m t ) + 2B(p 2 t ,m g ,m h ) 

+ 4£ (p? , m t , m fi ) + A^igt i C (p f 1 , -p t , m t , m~ g , m^ ) + 2 (m| - m| ) C (p f 1 ,-p t ,m t ,m g , m t ~ 2 ) 

F 3 F = —5—5 { 2m 2 (pf , m § , ) - B (p 2 , m^ , m^ )] + n gt i (jigti - 4m? ) C (p f , -p* ,m t ,m g , mf ) 

- (VgtiVfrz - 2m tVgti ~ %™%lign)C{pi v -Pt,m t ,m g ,mz 2 )} 
Fi = -2efj, gt iB(p 2 i ,m g ,m t ) + 4(m| - m 2 ^)B{p 2 g , \,m g ) 

+ 2m 2 t [B(p 2 t , m g , mz 2 ) - B(p 2 t , m g , m^)] + Am 2 g B(p 2 g , m t , m^) 

+ 4m? (m? - m\^C(p ix , -p t , m t ,m g , ) - 2m? (m^ + m| - 2m 2 t )C{p ii , -p t , m t , m g ,mz 2 
F£ = AeB{p 2 h ,m g ,m t ) - ±B{p 2 g ,\m g ) - AB{p 2 g ,m u m h ) 

- 4/Xg t iC(p t ~ i , -j?t, m u m g ,mi x ) - 2(m| - mf^C^ , -p t , m t , m g ,mz 2 ) 

F£ = -^^\ 2 rn 2 [B{p 2 ,m g , m~ t ) - B(p 2 ,m g ,mz )} + ^ti(4m 2 - n gtl )C(pz , -p t , mt,m 5 , m t ~ ) 

+ (HgtiHgt2 ~ 2m?/igii - 2m 2 fi gt2 )C(pz 1 ,-p t , m t , m g , m t - 2 )| 
F B = N c [/xtigC^ , -p t , m fl , A, m*) - figtiCip^ , -p t , m g ,m t ,X) - VigtCip^ , -p t , A, m ?1 , m § )] 

- 2C F /j, t igC(pz 1 , -Pum~ tl , A, m t ) 

corrections from real-gluon radiation: 



a., 



47r 2 m 



l-Msl 2 



+ 



*i 



(m fi - m t )/^ - mt/t lt - m g J g5 - 



47r 2 m + -, iV c 



^1 



+ 



a 2 C F 



irm 



C F lf + AT C /J 



re-normalization of the coupling constant: 



A/", 



n. 



47T 



- - 7£ + log(4vr) - log ( 



ll lr 2 A7 2 1 



finite shift of the Yukawa coupling relative to the gauge coupling in MS: 



AT, 



J\f< 



47T 



-N c - C F 
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decoupling of the heavy flavors from the running strong coupling: 



In these expressions the invariant integrals are defined as 

1 f d 3 p t d 3 P gd 3 q / \ (2gp A )(2gp B )--- 



MB ... 1 f d 6 p t d 6 pgd 6 q / 
A(m ) = /i 2e 



B(p 2 ,m a ,m b ) = /x 2e y 



(27r) n g 2 — m 2 



(2qp a )(2qp b ) 



(2tt) w [9 2 -m2][(g+p)2- m 2] 



C(pi,p 2 ,m a ,m b ,m c ) = fi 



2c 



d n q 



(27r) n [q 2 - m 2 a }[{q + p\) 2 - m 2 \[{q +pi + P2) 2 - m- 2 ] 
B(p 2 ,m a ,m b ) = dB(p 2 ,m a ,m b )/dp 2 

Note that in the angular integrals I^ b B "' from real-gluon radiation we do not add a minus sign 
for every term (2qpi), in contrast to Ref. |I"6|] f| In the numerical analyses we have inserted the 
mass of the decaying particle for the renormalization scale fiR. The scale parameter /i accounts 
for the correct dimension of the coupling in n dimensions, cancelled together with the 1/e poles. 
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Figure 1: The Feynman diagrams for the self-energies: (a) self-energy of the stop particles, 
including the mixing due to the second and third diagram; (b) top-quark self-energy; (c) 
gluino self-energy [including fermion- number violation]. 
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Figure 2: The dependence of 9(Q 2 ) on the renormalization scale Q. The normalized 
shift is shown relative to 9(m~^): [9{Q 2 ) — 9(m 2 ^)}/9(m 2 ^). The input mass values are the 
same as for the stop decay to gluinos: m.1/2 = 150 GeV, itlq = 800 GeV, Aq = 200 GeV, 
H > 0, for which the leading-order mixing angle is given by 1.24 rad. The minimum of 
the correction corresponds to the threshold Q = rrig + m t in the scalar integral. 
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Figure 3: (a) Born diagrams for stop and gluino decays; (b) vertex corrections for stop 
decays; (c) real-gluon emission for stop decays. The corrections to gluino decays can be 
obtained by rotating the diagrams. 
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Figure 4: The SUSY-QCD corrections to the decay of heavy stop particles into top quarks 
and gluinos: (a) the masses of the particles [in GeV] as a function of the common gaugino 
mass mi/2] (b) the decay widths in leading order (dashed curve) and next-to- leading order 
(solid curve). Input parameter set: m = 800 GeV, A = 200 GeV, \i > 0. [The kink 
at the threshold for g — > tii can be smoothed out by inserting the finite widths of the 
particles.] 
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Figure 5: The SUSY-QCD corrections to the decay of gluinos into light stop particles and 
top quarks: (a) the masses of the particles [in GeV] as a function of the common gaugino 
mass mi/ 2 ; (b) the decay widths in leading order (dashed curve) and next-to- leading order 
(solid curve). Input parameter set: m = 400 GeV, A = 200 GeV, ji > 0. 
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Figure 6: The SUSY-QCD corrections to the decay of light stop particles into top quarks 
and the possible neutralino eigenstates [see also Ref. ||]: (a) the masses of the particles 
[in GeV] as a function of the common gaugino mass m\ii\ (b) the decay widths in leading 
order (dashed curve) and next-to-leading order (solid curve). Input parameter set: m = 
50 GeV, A = 100 GeV, \i < 0. Only the decays into the two lightest neutralinos are 
kinematically allowed for this parameter set. 
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Figure 7: The SUSY-QCD corrections to the decay of light stop particles into bottom 
quarks and chargino eigenstates [see also Ref. ||]: (a) the masses of the particles [in 
GeV] as a function of the common gaugino mass m^; (b) the decay widths in leading 
order (dashed curve) and next-to-leading order (solid curve). Input parameter set: m = 
50 GeV, A = 100 GeV, fi < 0. 
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Figure 8: The SUSY-QCD corrections to the decay of top quarks into light t\ particles and 
light neutralinos [see also Ref.||]: (a) the masses of the particles [in GeV] as a function of 
the common gaugino mass m^; (b) the decay widths in leading order (dashed curve) and 
next-to-leading order (solid curve). Input parameter set: m = 250 GeV, A = 800 GeV, 
fi > 0. Only the decays into the two lightest neutralinos are kinematically allowed for 
this parameter set. 
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